Crossover from classical to random-field critical exponents in As-doped TDVO4 
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Using birefringence techniques we have measured the critical exponents /3, 7, and 8 in As-doped 
TbVCU, a structural realization of the random- field Ising model where random strain fields are 
introduced by V-As size mismatch. For pure TbVCU we observe the expected classical critical 
exponents, while for a mixed sample with 15% As concentration our results are /3 = 0.31 ± 0.03, 
7 = 1.22±0.07 and 8 — 4.2±0.7. These values are consistent with the critical exponents for the short 
range pure Ising model in three dimensions in agreement with a prediction by Toh. The susceptibility 
data showed a crossover with temperature from classical to random field critical behaviour. 



Determination of the. critical properties of the ran- 
dom field Ising modelErtJ (RFIM) has shown encouraging 
progress recently, both experimentally and theoretically, 
after many years of uncertainties. Measurements on di- 
lute antiferromagnets in a field, the most frequently stud- 
ied realization of the RFIM, have been difficult to obtain 
because dilution of the magnetic species inhibits equili- 
bration close to the critical temperature. However recent 
experiments by Slanic et ala on samples with only 7% 
dilution have shown equilibrium behaviour through the 
transition temperature, thus allowing confident determi- 
nation of the critical exponents. In the last few years, 
theoretical investigations making use of a variety of an- 
alytical and computational techniques have led to pre- 
dictions of the critical exponents of the RFIM that show 
reasonable consistency with, each other. A recent analysis 
by Fortin and HoldswortbO supports earlier suggestions 
that the critical exponents for the RFIM in three dimen- 
sions (d = 3) are those of the Ising model for reduced 
dimension d! = 1.5. Nevertheless dimensional reduction 
has not been rigorously proved, and even if the effective 
dimension d! — 1.5 is correct the uncertainty in some of 
the calculated critical exponents is quite large. Likewise 
the experimental situation is still not satisfactory, since 
the critical exponents measured by Slanic et ala have 
substantial uncertainties and are only partially consis- 
tent with theory. In addition, measurements of specific 
heat critical exponcntsEm in dilute antiferromagnets also 
appear to disagree with theory. 

For another realization of the RFIM, where random 
strain fields are generated by substitutional impurities in 
crystals undergoing structural Ising (Jahn- Teller) tran- 
sitions, the results to date have also not been conclu- 
sive. Random fields due to As/V substitutions in DyVC>4, 
which has d = 3 Ising exponents, appeared to increase the 
susceptibility exponent 7 as expected but had no effect 
on the order parameter exponent f3.u The interpretation 
of the effects of random fields in the As-doped DyVC>4 
system is complicated by the fact that the true critical 
behaviour of pure DyVC>4 should be classical due to the 
long range strain coupling. However because of the rel- 
ative weakness of the long range to short range interac- 
tions, classical exponents are not observable at accessi- 
ble temperatures \t\ > 1(T 2 , where t = (T - T D )/T D 



is the reduced temperature,™ leading to uncertainty on 
what the effects of the random fields would be. We 
have therefore extended these experiments to the related 
TbV04/TbAsC>4 system where the critical bekaTpsur of 
the pure compounds is unequivocally classicalMHla, and 
searched for changes in critical behaviour in mixed crys- 
tals due to random fields. Since this system starts from 
a different universality class, the results cannot be com- 
pared directly with results from dilute antiferromagnets, 
but it is an important system that can independently 
test theoretical models and predictions of random field 
effects. In contrast with the large number of theoretical 
investigations of the short-range RFIM, the literature on 
the random-strain version of the RFIM-is very limited, 
consisting primarily of a paper by Toh.113 In this paper, 
Toh compares the random-strain RFIM with long range 
forces to that of the short-range RFIM under a renor- 
malization group analysis. The main result is that the 
critical exponents should change from classical values to 
values that are close to those of the d — 3 pure short- 
range Ising model. 

The lowest 4/ electronic levels in TbVC>4 consists of 
2 singlet states ~ 18 cm -1 apart and a non-Kramers' 
doublet approximately halfway between. Coupling be- 
tween the doublets and lattice distortions leads to a Tb 
ion-ion interaction of the Ising formlij and a tetragonal- 
orthorhombic phase transition at temperature Td ■ Since 
the Tb ions are coupled predominantly to k ~ acoustic 
phonons and to bulk strains, the ion-ion interaction is 
very long range. The order parameter is the macroscopic 
strain a — b where a and b are basal plane unit cell pa- 
rameters in the orthorhombic phase. The orthorhombic 
distortion gives rise to birefringence, An, which is pro- 
portional to a — b (at least to a good approximationE£l) . 

The full Hamiltonian for the coupled clcctron-phonon 
TbVC>4 system in an external magnetic field can be writ- 
ten as 



H = -I E J «<°i - ^Ed + r « )°? - B E 



(i) 



f)rf and a z , a x , r z and t x are 
Jn describes the ion-ion inter- 



where mf = jg/ZBQ 
Pauli type operators^., . , , 
actions, 2e is the high temperature splitting between the 
outer singlets and B is the magnetic field applied along 
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the x (or a) axis (i.e. along the 110 direction). The field 
B is able to induce an orthorhombic distortion because 
of the strongly anisotropic Tb magnetic moment in the 
orthorhombic phase. In the mean field approximation, a 
Landau expansion of the free energy, with_£ = 0, shows 
that B 2 /T is effectively an ordering fieldJIZI and this is 
supported by the experimental data that follows. For B 
and e small, TbV04 is well described by an Ising model 
Hamiltonian. As the mode softening at the transition in 
this type of system is anisotropic!!!!, classical critical ex- 
ponents are expected, and observed, rather than d = 3 
Ising exponents. 

In the mixed compound, Tb(As 2; Vi_ a ;)04, a fraction x 
of the V atoms are replaced by As atoms, generating ran- 
dom, static strain fields, one component of which has the 
right symmetry to couple to the order parameter xB For 
e = B = the Hamiltonian has the form of the RFIM, 

ij i 

where hi is the random local strain field which is expected 
to have a Gaussian distribution about h = 0. In his anal- 
ysis of this type of random field system with anistropic 
mode softening, Toh predicts changes to the upper crit- 
ical dimension, thus modifying the values of critical ex- 
ponents at d = 3. 

Crystals of Tb(Asj;Vi_ x )04 with impurity concentra- 
tions of x — and 0.15 were prepared using the flux 
growth method at the University of Oxford. The crys- 
tals were cut and polished perpendicular to the c axis, 
with thickness of about 1 mm. They were mounted in a 
strain-free manner with the c axis horizontal in a helium 
optical cryostat. The crystal could be rotated about a 
vertical axis, allowing alignment of the c axis parallel to 
the laser beam and one of the a axes parallel to a hori- 
zontal magnetic field. A circular aperture in the sample 
holder of about 1 mm in diameter limited the sampled 
area to a small region of the crystal, and thus reduced 
effects due to any inhomogeneous composition, tempera- 
ture and ordering field that may be present in the crystal. 

The light source for birefringence experiments was a 
HeNe laser operating at 543.5 nm, a wavelength thai 
should give reasonable birefringence in this crystal.Ell 
Photoelastic modulation and lock-in detection allowed 
sensitive measurements of the phase shift <p due to the 
orthorhombic distortion. Adjustment of a Babinet com- 
pensator ensured that the detector output was propor- 
tional to <f>. In a typical experiment, the data acquisi- 
tion system brought the sample to each desired tempera- 
ture, waited up to 15 minutes for equilibration, and then 
ramped the magnetic field up and down while recording 
field, light intensity, and temperature. 

Below Td, the orthorhombic distortions are equally 
likely to be in the 110 or 110 direction, leading to the 
formation of twinned orthorhombic domains separated by 
domain walls. To avoid the cancellation of birefringence 
due to multiple domains, the crystal is forced into a sin- 
gle domain by applying a magnetic field that favours the 



distortion axis parallel to the magnetic field.c3c3 Thus 
to determine the birefringence, and hence /3, below Td, 
we recorded data over a range of magnetic fields and ex- 
trapolated the high magnetic field data back to zero field. 
With the rather small magnetic fields available (< 0.35 
T) it was sometimes difficult to achieve a single domain 
at temperatures close to Td , at |T - Td| < 0.2 K for 
x = 0.15, but this became progressively easier further 
away from Td- For the x = sample where pinning 
is presumably weaker, single domain structure at |T - 
Td| < 0.75 K was relatively easy to achieve. For this rea- 
son, the critical isotherm exponent, <5, obtained from the 
dependence of the induced birefringence on ordering field 
at the transition temperature, is more reliable for pure 
TbVC>4 than that for the mixed sample. Above Td, the 
change in induced birefringence resulting from a change 
in the ordering field gives the susceptibility and hence 
7. Although the exponents for pure TbV04 are known 
to be classical,E£l, their measurement provides a compari- 
son with the mixed sample where modified critical expo- 
nents are expected. In experiments, Td appeared to vary 
slightly from run to run because of effects such as mount- 
ing strains and temperature gradients. To minimize the 
effects of a variable Td on the results, efforts were made 
to measure all exponents of one particular sample on the 
same run. In the power-law fits to the data, Td was 
chosen to optimize both the (3 and 7 fits simultaneously. 
The consistency in the results of repeated experiments 
and the quality of the fits give confidence in the results. 

The birefringence with increasing and decreasing or- 
dering fields for various temperatures below Td in the 
x = 0.15 sample is shown in Fig. |l|. Hysteresis observed 
in the low field region is attributed to pinning of the mul- 
tidomain structure. The data from the higher field region 
where hysteresis is absent were extrapolated to determine 
the zero field birefringence. A log-log plot of both the 
x = 0.15 and x = data, fitted to a power law of the form 
An tx \t\P is shown in Fig. |. Using T D = 29.26 ± 0.03 K 
and data in the range 0.005 < \t\ < 0.027, we obtained 
a value of/3 = 0.31 ±0.03 for the mixed sample. At 
larger \t\ we found no convincing evidence for crossover 
behaviour towards the classical exponent. For the pure 
sample, we obtained /3 = 0.46±0.06 for T D = 32.32±0.04 
K from data in the range 0.004 < |i| < 0.037. 

Figure || presents susceptibility data for selected tem- 
peratures above Td- At temperatures close to Td, these 
slopes deviate from linearity, attributable in part to the 
presence of non-zero birefringence at zero field, probably 
caused by internal strainalj. The susceptibility exponent 
7 was found by fitting the various values of x(T) to the 
relation x~ l ^ ■ Figure || shows a log-log plot of the 
susceptibility \ versus the reduced temperature \t\ for 
x = and 0.15 samples. As can be seen from Fig. |J, 
the log-log plot for x = 0.15 shows two linear fits with 
the data closer to Td having a larger slope than that for 
data further away. Data in the range 0.027 > \t\ > 0.01 
(closest to Td) were optimized first to a linear fit giving 
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T D = 29.26 ± 0.03 K and 7 = 1.22 ± 0.07. A linear fit to 
data in the range 0.100 > \t\ > 0.021 (further from T D ) 
with Td unchanged yields a value of 7 = 0.89 ± 0.03. We 
did not attempt to fit these data to a crossover function, 
but they suggest a crossover temperature near \t\ — 0.024 
(29.96 K). The log- log plot of the susceptibility versus 
reduced temperature for pure TbV04 in Fig. || does not 
show a crossover effect. A power law fit in the range 
0.058 > \t\ > 0.005 gives T D = 32.32 ± 0.04 K and 
7 = 0.92 ±0.07. 

Nonlinearity in the dependence of the order parameter 
on the ordering field at temperatures below Td can be 
attributed to the pijogressive detwinning of the sample 
with increasing fieldc-a (i. e. changing from a multido- 
main to a single domain structure). The critical isotherm 
exponents d were extracted from the field-induced bire- 
fringence data at temperatures closest to that of the pre- 
viously determined transition temperatures for x = 0.15 
and x = samples. Only the higher field An data were 
fitted to the power law (B 2 /Td) 1/<5 . For the x = 0.15 
sample, we obtained a value of S = 4.2 ± 0.7 at T = 
29.24 ± 0.03 K and for the pure sample 5 = 2.6 ± 0.4 at 
T = 32.33 ± 0.03 K. The uncertainty in locating T D is 
incorporated in the uncertainty in 6. 

Our values of the critical exponents for pure TbV04 
agree satisfactorily with the values f3 = 0.5,7 = lj an d 
(5 = 3 expected for a mean field system. For the random- 
field sample our values /3 = 0.31 ± 0.03, 7 = 1.22 ± 0.07, 
and 5 = 4.2 ± 0.7 are in good agreement with the expo- 
nents for the d — 3 Ising model (3 = 0.33, 7 = 1.24, 
and 5 = 4.8. n 

Toh's predictionsli3 for the effects of random fields on 
the exponents in this type of system are therefore well 
supported by our results. Toh noted that while the 
anisotropic strain interactions in TbV04 reduce the up- 
per critical dimension d* from 4 for the standard Ising 
model to 2, resulting in classical exponentsEJ, the intro- 
duction of random strain fields raises d* from 2 to 4 again, 
leaving the critical exponents the same as for the stan- 
dard Ising model in d = 3. In view of the focus in recent 
theoretical analysis on dimensional reduction it is of in- 
terest to examine Toh's analysis and our results in terms 
of changes in effective dimension d 1 instead of in d* . Thus 
for the Ising model we regard the upper critical dimension 
to be fixed at 4; in pure TbV04, the anisotropic strain in- 
teractions presumably raise the effective dimension from 
3 to 5, giving classical exponents as before. If changes in 
critical properties due to random fields can be described 
by dimensional reduction, the consensus predictiontffl is 
that random fields reduce the effective dimensionality by 
2 — 7], where rj is the exponent for the decay of magneti- 
zation cotselations. In the present case, since r\ = 0.032 
for d = 3,c3 the dimensionality should be reduced by 1.97 
or essentially 2, resulting in d' = 3. Hence our results, 
within their accuracy, are consistent with dimensionality 



reduction by random fields of 2 — r\ ~ 2 for this system. 
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FIG. 1. Changes in phase <p (birefringence) in x — 0.15 sample at selected temperatures below Td for increasing (lower 
points) and decreasing (upper points) ordering fields. The lines join data points. 
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FIG. 3. Changes in birefringence versus ordering field in x — 0.15 sample at selected temperatures above To- Slopes in the 
small-field limit give susceptibilities. 



G 




-2.4 -2.0 -1.6 -1.2 

Logt 



FIG. 4. Log-log plots of susceptibility versus reduced temperature for x — and 0.15. The data for the mixed sample show 
a crossover from classical to random-field critical behaviour. 
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